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Abstract 

In this review we discuss novel phases (quantum critical points occurring only 
for T = 0) in (2+l)-dimensional Abelian gauge theories, which may serve as proto- 
types for studying the physics of nodal excitations in (underdoped) high temperature 
superconductors. We pay particular attention to describing the existence of novel 
phases related to the anomalous breaking of certain symmetries, including uncon- 
ventional superconducting properties. Although some of the results are rather old, 
we present them here from a physical perspective not discussed previously in the 
literature. In this respect we also explore the (in) applicability of some non pertur- 
bative approaches on the existence of a critical number of fermion flavours for chiral 
symmetry breaking in our specific models, which is an issue closely related to the 
existence of the novel phases mentioned above. 



1 Introduction 



Relativistic gauge field theories in (2+l)-dimensioris play an important role in at- 
tempts to understand the dynamics and the associated phase diagrams of planar 
doped antiferromagnets 1J, and through this to shed light on the still elusive micro- 
scopic theory of high temperature superconductivity EJ IH El E] ■ The relativistic 
nature of the pertinent excitations arises when considering the effective (continuum) 
field theory that governs the degrees of freedom near the nodes of the Fermi surface 
or points where the superconducting gaps vanish. Indeed, high temperature super- 
conductors are experimentally known to be strongly type II d-wave superconduc- 
tors, which have their superconducting gaps vanishing at some points in momentum 
space. One popular approach to the physics of high T c is to linearise a spin-charge 
separating theory [7] of holon-spinons around such nodes, resulting in relativistic 
gauge theories of fermions and bosons coupled to gauge fields. This type of theories 
provide a description of the effective interactions in the ground state between the 
fundamental degrees of freedom; in this picture the physical electrons are not con- 
sidered as fundamental particles, being instead composites (bound-states) of spinons 
and holons [7j. 

There are two main approaches to this problem, distinguished by the spin and 
statistics properties of these electron constituents: (i) the slave-fermion approach [HI 
12], in which the holons are viewed as Dirac fermions, electrically charged, and the 
spinons as bosonic neutral fields, and (ii) the slave-boson approach jl], according 
to which the spinons are viewed as neutral fermions and the holons as charged 
bosons. The two approaches are supposed to be physically equivalent, by means of 
appropriate non-Abelian bosonisation in a path integral continuous formalism [U]. 
However, within each approach there are various models considered in the literature, 
which may not be physically equivalent. The associated continuous effective theories 
depend crucially on the way the continuum limit is taken. This is an important point, 
often ignored in the recent literature on the subject. The upshot of the present work 
is to critically examine such models from the point of view of recent non-perturbative 
results on the symmetry structure of three dimensional gauge theories, and argue on 
the existence of novel quantum critical phases; the latter, although known from the 
early literature on three-dimensional gauge theories, seem not to have been taken 
into account in the recent (condensed-matter) studies on the subject. 

In this work we shall not give a detailed description of condensed matter results, 
but rather concentrate on the comparison of the various effective continuum gauge 
models which have been argued to describe the dynamics of nodal excitations in 
doped antiferromagnetic materials. Even though the details of the microscopic the- 
ory are important for arriving at specific models they will not be crucial for our 
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analysis. We refer the interested reader to the relevant literature for details on these 
issues. For our purposes we shall simply associate the dynamics of the spin-charge 
separation in condensed matter systems to effective continuum field theories, which 
are variants of three dimensional quantum electrodynamics (QED3) with 2N flavours 
of two-component Dirac fermions. To that end we shall adopt the simplest version 
of the spin-charge separation leading to an Abelian statistical gauge interaction 
between spinons and holons. 

In what follows we will discuss the various models, with emphasis on those phys- 
ical properties which indicate the appearance of novel phases, not discussed in the 
recent literature. In addition, we will critically examine the possibility of restrict- 
ing or excluding such phases by resorting to recent non-perturbative arguments ^U] 
related to the chiral symmetry structure of (three-dimensional) gauge theories. In 
particular, we shall argue that, due to a variety of reasons, the constraints of JUj 
are not applicable to the models considered in this work. 

The outline of the article is as follows: In section 2 we discuss the effective 
field theories of relevance for the underdoped phase of planar antiferromagnets. We 
put emphasis on the unconventional phases, occurring either at T = or at most 
up to very low (mK) critical temperatures, which may be superconducting in the 
(Landau) sense of having massless poles in the electric current-current correlator. 
Such phases are characteristic of non compact theories, and are absent for compact 
ones, where no massless poles appear in the aforementioned correlator. We discuss 
such phases for all variants of three-dimensional Abelian gauge theories of interest 
to condensed matter. In section 3 we discuss critically some non-perturbative con- 
jectures of [10. on the chiral symmetry breaking structure of gauge theories. If this 
type of conjectures (which are shown to hold in most four-dimensional examples) 
were valid, this would challenge the existence of quantum-critical phases for some 
of the models, but not for all (specifically the r 3 -QED model of [2] escapes this con- 
straint). However, we present a variety of reasons as to why the inequality of [TP] 
may not be applicable to the class of models that we consider. These reasons range 
from infrared (IR) infinities, which invalidate some of the counting arguments, to 
potential misconceptions related to the ultraviolet (UV) behavior of QED3 and in 
particular its (non)-display of asymptotic freedom. This last notion, as well as the 
issue of the existence of a critical number of fermion flavours for a breaking of chiral 
symmetry are examined from the point of view of Schwinger-Dyson (SD) equations 
in the context of a specific approach in section 4. In particular, the appropriate 
definition of what we call a (dimensionless) effective charge in QED3 is given, and a 
non-linear SD equation is derived for the semi-amputated vertex, which constitutes 
the natural non-perturbative generalization of the aforementioned concept. The 
resulting integral (SD) equation governs the dynamical evolution of the coupling, 
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and can be solved by transforming it to an equivalent differential equation of the 
Emden-Fowler type. Its solution indicates that unless some type of mass is dynam- 
ically generated the coupling displays unbounded oscillatory behaviour in the IR. 
When the coupling obtained from the semi-amputated vertex is inserted into a stan- 
dard gap equation it is found that chiral symmetry does take place and a fermion 
mass is indeed generated dynamically A subsequent comparison of our results to 
those obtained using large-A techniques is made. Our main conclusion from this 
SD analysis is that, although there exist regions of allowed values of this effective 
charge for chiral symmetry breaking to occur, nevertheless there is no restriction 
on the allowed number of fermion flavours. Conclusions and outlook are presented 
in section 5. Finally in an Appendix we review briefly some of the salient facts of 
the pinch technique (PT) • The latter is a systematic method for construct- 

ing gauge-independent off-shell Green's functions, and could in principle lead to a 
manifestly gauge-invariant truncation of the SD series. 

2 Novel effects in nodal liquids 

In this section we present arguments pointing towards the existence of anomalous 
superconductivity and new quantum critical points in nodal liquids at T = 0. 

2.1 An effective field theory approach 

Relativistic superconductors in (2+l)-dimensions, where superconducting properties 
are due to some anomalous graphs in the effective continuum gauge theory, have 
been first considered in [2], as a way of explaining parity- invariant high-T c planar 
superconductors with either nodal points in their Fermi surface, or Fermi surfaces 
consisting of small isolated spherical pockets (or points). The basic idea behind such 
models is spin-charge separation [7j, according to which the fundamental degrees of 
freedom in the ground state of doped antiferromagnets (AF) are not the ordinary 
electrons, but two kinds of confined constituents, one with possibly fractional spin, 
but electrically neutral, termed spinon, and the other with electric charge, and 
possibly spinless, termed holon. 

The interest in nodal liquids has been recently revived in view of experimental 
results confirming the o?-wave nature of high temperature superconductors, as well 
as the fact that the Fermi surface of these materials in the so called underdoped 
phase consists of four nodal points. In the recent theoretical literature a lot of works 
have appeared linearising the spin-charge separating theory about such nodes, in an 
attempt to understand the dynamics of the underdoped cuprates. 
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There are in general two formal approaches to describe the field theory of spinon 
and ho Ions: (i) the so-called slave-boson approach 0] treats the spinons as neutral 
fermions and the holons as charged bosons, and (ii) the so-called slave-fermion 
approach [HI2]) which treats the spinons as neutral bosons and the holons as charged 
fermions. It turns out that the field-theoretic description of the dynamics around a 
nodal point in the Fermi surface is best accomplished through an effective relativistic 
gauge field theory, where the gauge field expresses statistical effective interactions 
between spinons and holons (essentially spin-spin interactions, due to the magnetic 
origin of superconductivity in the spin-charge separation approach [7j). 
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Combining Nodus 1+3 yield 
one four component 
spinor, whilst combining 
nodes 2+4 yields another 



Figure 1: In the QED 3 models of J<3 [3J/ the four component spinors (spinons) of the 
continuum effective theory are constructed by combining nodes along the diagonal; 
their flavour index expresses the existence of two such pairs. In this scheme there are 
variations as to the precise nature of the components of the four- component spinors, 
which may lead to important physical differences JijJJ/. 



The various models existing in the literature may be classified in two major 
categories, depending on the nature of their fundamental excitations: 

• (A) In the slave-boson approach of E], the effective nodal theory of the 
spinon part is nothing but QED 3 , coupled to charged boson degrees of free- 
dom, corresponding to the holons. In this model the gauge interactions are 
assumed to be Abelian, non-compact interactions. The way the four compo- 
nent continuous spinors, representing the spinons, are constructed is by means 
of combining nodes as shown in figure ^ There are variations on the precise 
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construction of the effective microscopic degrees of freedom, which may lead 
to important physical differences of the resulting continuum theory, especially 
as far as the insulating phase of the materials is concerned We shall not 
dwell upon these differences here, given that they have little or no bearing on 
the issue of interest, namely the appearance of novel quantum critical points; 
the latter seem to be common to all approaches mentioned above, and will be 
the focal point of the discussion that will follow. 
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Figure 2: In the r 3 — QED model of Jjfl, the four component spinors of the contin- 
uum theory (holons ) are constructed by combining lattice points, and then there is a 
"colour" r 3 index coming from the Antiferromagnetic sublattice structure. 



• (B) In the slave-fermion approach considered in |2] the four-component spinors 
of the resulting effective nodal theory represent the electrically charged holons, 
while the spinons are CP 1 bosons (magnons). The four component holons are 
constructed by combining appropriately the antiferromagnetic sublattice struc- 
ture, as indicated in figure 121 The resulting continuum theory is the so-called 
r 3 QED 3 model, in which the statistical gauge field couples with opposite cou- 
pling to the two 'colours' of the holons, representing ho Ion excitations in each 
sublattice. The gauge interactions of the model may be actually embedded in 
the full non-Abelian SU(2) gauge group structure [T3|, which results in the pos- 
sibility of having compact Abelian gauge groups . As we shall discuss below, 
there are important physical differences, between compact and non-compact 
cases. 

The two approaches mentioned above are not physically equivalent, due to the 
different ways by which the four component spinors are constructed, but also due 

1 For an alternative SU(2) formulation of the doped t — J model, within the slave-boson frame- 
work, see |14| . 
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to the distinct origin of the statistical gauge fields. One should not confuse this 
issue with the one in which the slave-boson and slave-fermion Ansatze are related 
by means of non-Abelian bosonization in a path integral approach . Each one of 
the models in (A) and (B) above can be bosonised or fermionised accordingly, but 
the resulting models will not connect (A) to (B). They will simply be the bosonised 
(or fermionised) versions of (A) or (B). This can be easily understood from the 
fact that, in the approach (B), for instance, the spinons and holons couple to the 
statistical gauge field with opposite couplings, due to the AF sublattice structure, 
a feature which will persist upon bosonization. This would imply simply that the 
[spinor spinons] of the fermionised approach derived as in jHj will be different from 
those of [UJ Ej, whose coupling with the statistical gauge field does not have the 
above-mentioned AF "color" structure. 

Despite this inequivalence, however, both approaches are still characterised by 
anomalously-broken symmetries, leading to quantum critical points at T = 0, al- 
though the broken symmetries are different in each case. The purpose of the fol- 
lowing two subsections is precisely to analyze these phases in some detail in the 
context of (A) and (B). Even though the possible existence of such phases has been 
advocated for quite some time [2] , in our opinion their relevance has not been 
fully recognized by the condensed-matter community. 

We shall argue in this work that their existence may have important implications 
for the entire phase diagram of the high-T c , and more general the doped AF. We 
commence our discussion from the case (B), which historically was the first (2+1)- 
dimensional relativistic model of gauge field theories to be used in the context of 
planar doped AF, generalising a (l+l)-dimensional spin chain model of [T3] . 

2.2 T3-QED in the slave-fermion spin-charge separation 

We start from a brief description of the results of [2] , which utilizes the slave-fermion 
approach to spin-charge separation. In this case, the electron operators are written 
as 2 : 

Ci,a = ^PiZ a , a = 1,2 (2.1) 

where ipi (i is a Lattice site index) denotes the spinless holons, which are Grassmann 
number on the lattice carrying electric charge but no spin, and z a are CP 1 magnons, 
representing the spinons, with a = 1, 2 a spin SU(2) index. There is of course the 

2 For the most part of this work we restrict ourselves to Abelian spin-charge separation Ansatze 
for brevity; one can carry over the discussion to non abelian cases, where the statistical gauge 
group is non abelian we shall make some comments on the physics of the non-abelian case 
later on in the article. 
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constraint of having at most one electron per lattice site 



2 



4a + 



ZiZi = 1 , (no sum over i) 



(2.2) 



a=l 



We first give a brief but comprehensive description of the continuum construc- 
tion presented in [2], where the interested may find more details. In that work, 
the antiferromagnetic nature of the underlying lattice defines two "colours" of two- 
component fermions; one obtains two additional flavours because of the usual lattice 
doubling (c.f. figure EJ). In terms of nodes these are the degrees of freedom around 
each node in the construction of In contrast, in the constructions of |HE!, which 
can also be considered in the context of T3-QED JE] one combines the two nodal 
points along the diagonal of the appropriate Fermi graph, to construct one species 
of four component spinor, and then one duplicates the species by taking into ac- 
count the combination of the other two nodes, as indicated in figure ^ However, 
in that case one should couple the two so-obtained four- component spinors to the 
gauge fields with opposite statistical couplings in order to get the r 3 -QED model, 
which allows the statistical gauge interactions to be viewed as forming a subgroup 
of the gauged spin SU(2). This is the important physical difference compared to the 
construction of jHl IS] , where one does not distinguish the sublattices, as far as the 
statistical gauge interactions are concerned. In that case, the latter are not directly 
related (embedded) to the gauged spin SU(2) group. 

The effective continuum relativistic nodal (2+l)-dimensional lagrangian reads |2j: 



where z are spinons (bosons), and \l/ c , c = 1, 2 are four-component fermions, with a 
colour index associated with the antiferromagnetic nature of the underlying micro- 
scopic lattice t — J model j2] . The gauge field a M is an abelian field expressing the 
effective spin-spin interactions among the fundamental excitations of the model, and 
is different from the real (external) electromagnetic field A^, to which the electrically 
charged holons \l/ c couple. Notice that, as a result of the antiferromagnetic nature, 

there is a r 3 = ( J J j coupling for the a-field, but an ordinary QED coupling 

for the A-field. Moreover, the effective "speed of light" of the relativistic continuum 
theories is the Fermi velocity of the node, which here is taken to be one. It is for this 
reason that the real electromagnetic speed of light , c has been kept explicit in the 
coupling in (|2.Hjl . In realistic models it is estimated that c = 10 4 in units of the 
fermi velocity of the nodes. This should be kept in mind throughout this work. The 
. . . in ()2.3|) express contact interactions among the fermions. From them, the most 



£ = 7 Tr|(^ - ia^z\ 2 + #(*7 M ^ + 7^3 + VAJ* + ■ ■ ■ 



(2.3) 
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relevant ones are four-holon interactions, which in (2+l)-dimensions are known to 
be relevant in a Renormalization-group sense ^7] . We shall come back to their role 
in connection with the claims made in [TO] in this context in the next section. 

We now remark that the bare mass of the holons \l/ c is zero, but a parity conserv- 
ing mass can be generated dynamically by means of the statistical gauge interactions 
a n P"%1I2]. The parity conserving mass is energetically preferred due to a theorem in 
[TU] for vector-like theories. In the absence of external fields A^, and in the massive 
phase of the spinons z, which can thus be integrated out producing kinetic Maxwell 
terms H(0)f£ v (Jl(k) is the one loop vacuum polarization), such a mass is responsible 
for a spontaneous breaking (i.e. by the ground state of the system) of the global 
fermion-number symmetry, generated by the current J fl = J^ c - The breaking is 
due to the anomalous graph of fig. El as discussed in detail in The T3 coupling 
is crucial to this effect. 

The resulting matrix element is: 

, , M p p . . 

< >~ M <W-pj (2-4) 



where M is the mass of the fermion The result (|2.4jl is perturbatively exact. 



Figure 3: The breaking of fermion number symmetry upon mass generation of 
fermions in the r 3 QED 3 model of L 2l. The figure denotes the matrix element of 
the fermion number current between the vacuum and one 'photon 7 state. The dark 
blobs denote loop corrections and the crossed-blob denotes an insertion of the fermion 
number current. 

Upon coupling to external electromagnetic fields A^, the anomalous matrix ele- 
ment (|2.4jl (c.f. fig. EJ) leads to superconductivity in the case of non compact gauge 
fields result of the existence of a massless pole in the electric current-current 

correlator (c.f. fig. HJ: 

< o| j„(ao J„(-*)|o >= ^ - O (1 + n^r 1 ( 2 - 5 ) 

where H(k) is the one-loop vacuum polarization graph. 
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Figure 4: In non compact gauge theories the anomalous graph of fig. d gives rise 
to massless poles in the current- current correlators, as a result of the exchange of 
the massless statistical gauge boson a^; upon coupling the system to an external 
electromagnetic potential, this results in superconductivity (Landau criterion). 

Strictly speaking this superconducting behaviour pertains only to zero temper- 
atures T = 0, since at finite temperature the masslessness of the statistical gauge 
field disappears, as a result of a plasmon mass in the longitudinal component 
a . In view of the above considerations therefore, the novel T = phase predicted 
here would result in a modification of the temperature-doping phase diagram of 
high-T c , by a quantum critical line at T = 0, as shown in figure should the r 3 - 
model describe correctly the underdoped cuprate phase. It has been argued, though, 
in [2] that despite the plasmon longitudinal mass, the screening of the magnetic field 
(Meissner effect), associated with the massless transverse component of the statisti- 
cal gauge field a^, is still valid at finite temperatures, up to a temperature in which 
the nodal holon mass gap disappears. Such a temperature has been estimated in 
|2U] to lie at most on the mK scale, i.e. much lower than the optimal- doping critical 
temperature of 100 K of high-temperature superconductors. 

2.3 QED3 in the slave-boson approach 

The issue we want to bring up in this subsection is that in the non compact gauge 
a-field corresponding T = anomalous behaviour characterises also the 

QED 3 models of jHl Ej, within the slave-boson approach, but the symmetry that 
breaks anomalously in this case is not the fermion-number symmetry, but a version 
of chiral symmetry (but not the standard chiral symmetry which breaks due to 
the parity-conserving mass generation for fermions). In this sense, the quantum 
critical point does not imply superconductivity, but an anomalous behaviour within 
an AF phase. This quantum-critical phenomenon has been missed in the relevant 
recent literature E], because the analysis is restricted only to T ^ 0. It is 
our conjecture, therefore, that in QED 3 theories of the pseudogap phase of high- 
temperature superconductors utilizing non compact gauge fields, there is always an 
anomalous behaviour in the spin sector at T = 0, which is realized via the loop 
graphs of figs. 13141 By the way, we call this behaviour "anomalous" since at tree 
level the result is zero, and appears only at one loop. 
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normal ? 
(non fermi liquid) 








Figure 5: The anomalous graph of fig. [3] results to a novel quantum critical phase, 
of Kosterlitz Thouless (KT) type (no local order parameter perturbatively) strictly at 
T = 0. Should the model of JJfl describe the underdoped cuprate phase, then, this may 
modify the phase diagram of high temperature superconductors, by connecting the 
optimal doping superconducting phase and the antiferromagnetic phases, as indicated 
by the dashed line. It must be stressed, though, that if the screening of the magnetic 
field lines still occurs for finite temperatures up to a critical T* much lower than the 
optimal T c ~ 100 K of high T c J^j, then the novel KT superconducting phase extends 
also inside the pseudogap phase in a region bounded by such small T*, as indicated 
by the dashed dotted line. 

Let us try to decipher this anomalous behaviour from studies of effective QED3 
models of the slave-boson approach to spin-charge separation developed in |H] • 
According to the latter, the electron operators are written as: 

Ci,a = Zi,aK a = 1,2 (2.6) 

where b are spinless bosons, carrying electric charge degrees of freedom, and rep- 
resent the holons, while £j iQ , are electrically neutral fermions, carrying spin, which 
represent the spinons, with a = 1, 2 a spin SU(2) index. There is of course the 
corresponding constraint (|2.2j) again, expressed in terms of the new variables in 

(EH). 

In jSJ E] an effective continuum field theory for £ has been constructed, which 
yields an interacting QED3 model of two four- component spinors \l/ constructed 
appropriately out of £ a , a = 1,2, interacting with a non compact statistical U{\) 
gauge field, expressing spin frustration responsible for nodes mixing. The lagrangian 
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reads: 



1 2 - 

Cqeds = —JIM + */(^ + + • • • (2.T) 

/=i 

where \&/ are electrically neutral four-component continuum spinors (spinons), and 
/ = 1, 2 runs over node pairs as in figure^ The massive \P phase is characterised first 

of all by a breaking of (global) chiral symmetries generated by 75 = i ( ^ ^ J and 
73 = % ( ^ J , in a four- component notation for spinors with an even number 



1 , 

of flavours, in a reducible representation of the Dirac algebra in (2 + l)-dimensions. 
Above, I denotes the 2x2 unit matrix. These are global symmetries, whose breaking 
will result in the appearance of massless Goldstone bosons, to be discussed in the 
next section. 

A dynamical mass generation (spin gap) for the spinons \I/ / has been argued 
in IB] to be related with properties of the pseudogap phase of high temperature 
superconductors. Most importantly, the massive spinon phase has been interpreted 
in [U] as implying that the antiferromagnetic (AF) phase is immediately succeeded 
by the superconducting phase, as indicated in figure El in the sense of a region of 
weak AF (spin density wave (SDW) phase) in which the critical temperature for the 
new phase is much smaller than the pseudogap critical temperature. This situation 
is opposite to the graph in figure of the model of [2], where the superconducting 
phase (for non compact gauge fields) enters the pseudogap region until the AF (Neel) 
state. The analysis of [B] did a low temperature analysis, attempting to generalize 
the conclusions down to the T = region. It is our purpose here to point out some 
important aspects of this region in QED 3 , not discussed in jB], but existing in the 
literature which like the T3-QED case, imply a discontinuity of the T = line 
(quantum critical line). 

Indeed, apart from the conventional chiral symmetry in theories with even num- 
ber of fermion flavours, discussed in jB], generated by the 4x4 matrix 73,75, there 
is another symmetry of ()2.7j) [3J 12], whose breaking is realized in the so-called 
Kosterlitz-Thouless mode, without a (perturbative) local order parameter, and hence 

Goldstone bosons. This symmetry is generated by the f 3 = i ( ^ ^ ^ The corre- 
sponding current generator is J^ 3 = Y^f=i ^ f^3l^ f an d, as shown in |3] its matrix 
element between the vacuum and the one photon state (as in fig. |3]) is non-zero in 
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KT SDW 



Figure 6: In the QED 3 model of [3J/ ; the anomalous graph of fig. [3J but for the 

chiral symmetry breaking current ^f^^ , results to a novel quantum critical phase, 
of KT type (no local order parameter perturbatively) which is a spin density wave 
(SDW) (weak AF) phase. This novel KT superconducting phase extends inside the 
pseudogap phase, and covers a region from the T = line to a critical temperature 
which is much smaller than the critical temperature separating the pseudogap from 
the normal phase. For any T ^ it becomes indistinguishable from the conventional 
chiral symmetry broken phase of 1$, but its T = anomalous behaviour, with a 
massless pole in the two-point correlator of the r^-chiral symmetry currents, pro- 
vides a distinctive discontinuous novel behaviour. Compare /contrast the situation 
with that in figure where one obtains an anomalous KT superconducting phase 
connecting the optimal doping superconducting (SC) phase with the AF phase. Here 
one obtains an anomalous SDW phase that connects the AF and (optimal doping) 
SC phases. 



the massive \I/ phase: 

, - , M p p , x 

< a„ W |0 >~ m ^j= (2.8) 

This also leads to Jl 3 current-current correlators with a massless pole for strictly a 
T = case, according to the corresponding graphs of fig. 0] (c.f. (I2.5jl ). Our claim 
is that such a KT T = phase, with a massless pole in the respective correlator of 
the generating currents defines a novel quantum critical point which corresponds to 
an unconventional AF phase in the model of |6,, or, in fact, any other QED3 model 
(within the slave-boson approach) of doped AF. As in the case of T3-QED, the critical 
temperature at which the chiral-symmetry-breaking spinon mass disappears, is also 
a critical temperature for this anomalous behaviour. This temperature is much 



12 



smaller than the critical temperature separating the pseudogap from the normal 
phases in the phase diagram. The distinctive feature of this anomalous symmetry 
from the standard chiral symmetry, which is used in the approach of 6\, is precisely 
its T = behaviour, where the current-current correlator for the spinon current 
^f 3 7 M \l/ has a massless pole that disappears at any finite T, due to plasmon mass. 
This feature defines the T = as a quantum critical line, corresponding to a KT 
breaking of the r 3 chiral symmetry, defined above: we may term this phase as KT 
SDW (c.f. figure EJ). For any T ^ the situation is like in j^j, but the T = 
quantum critical behaviour is distinct, and 'discontinuous' from the T^O case, and 
has been missed in the recent condensed matter literature [HIE]- 



2.4 Compact Gauge Theories 

In the above analysis the statistical gauge field had been assumed non compact. 
However, the spin-charge separation Ansatz might be extended to a non-Abelian 
form ^3] by exploiting appropriately an (approximate) particle-hole symmetry of 
the t — j Hamiltonian for underdoped cuprates: 

X*. = = { c , _ c » J ( = { _ 4 ^ ) i { Z2 - z J_ (2-9) 

where the fields z a ^ obey canonical bosonic commutation relations, and are associ- 
ated with the spin degrees of freedom (magnons), whilst the fields ip a ,i, a = 1, 2 have 
fermionic statistics, and are assumed to create holes at the site i with spin index a 
(holons). The Ansatz ()2.9|) has spin-electric-charge separation, since only the fields 
ip carry electric charge. It is a slave-fermion Ansatz since the holons are fermionic. 

The Ansatz has a local SU(2) symmetry, if one defines the transformation prop- 
erties of the z fields to be given by left multiplication with the SU(2) matrices, and 

*/3 



those of the ifj^g matrices by the left multiplication 



£a/3 - h&f, (2.10) 

In this representation, the gauge group SU(2) is generated by the 2x2 Pauli ma- 
trices. 

The Ansatz (J2.9|) possesses an additional local Us(l) 'statistical' phase symmetry, 
which allows fractional statistics of the spin and charge excitations. This is an ex- 
clusive feature of the three dimensional geometry. This is similar in spirit, although 
implemented in an admittedly less rigorous way, to the bosonization technique of 
the spin-charge separation Ansatz of ref. jU], and allows the alternative possibility 
of representing the holes as slave bosons and the spin excitations as fermions. 
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In addition, as a consequence of the fact that the fermions ip carry electric charge, 
one has an extra U em (l) symmetry for the problem. 

To recapitulate, the above analysis, based on the spin-charge separation Ansatz 
(12. 9 j) which allows spin flip, leads to the following local-phase (gauge) group structure 
for the doped large-C/ Hubbard model: 

G = SU(2)®U s {l)®U em {l) (2.11) 

where the second U em (l) factor refers to electromagnetic symmetry due to the elec- 
tric charge of the holes. This symmetry appears as a hidden symmetry of the effective 
holon and spinon degrees of freedom obeying the Ansatz ()2.9|) . 

The statistics-changing Us(l) group is strongly coupled, given that in the re- 
sulting effective lagrangian appears without a Maxwell Kinetic term, and hence it 
may be considered of (formally) infinitely string coupling g s — > oo. In practice, the 
coupling is cutoff by an effective cutoff which may be defined by the highest scale 
in the problem (e.g. the strong Hubbard coupling U in Hubbard models with a 
strong repulsion U ^> 1 so as to implement the one-electron per site constraint, or 
the Heisenberg interaction Jint — J models [13J). Such strong Abelian groups may 
be responsible for dynamical mass generation of fermions (holons). 

The effective continuum theory obtained from (|2.9|) around the nodes of the 
underdoped cuprates consists of Dirac spinors, representing the holons, as well as 
z magnons, interacting via the non-Abelian gauge interactions. The holons are 
coupled to the full Non-Abelian group G, while the z magnons, being electrically 
neutral, couple only to the SU(2) ® Us(l) subgroup. 

The dynamically generated mass for the holons is parity conserving as a result of 
the vector- like nature of the interactions in the simplest model considered [T9~] . As 
discussed in detail in [T3], such parity-conserving mass term break dynamically the 
SU(2) subgroup, due to the fact that the parity conserving mass term for fermions 
transform like triplets under this SU(2). 

As a result of the breaking SU(2) — > U c (l), two of the gauge bosons of SU(2) 
acquire heavy masses, and decouple from the effective low-energy theory. One is 
then left with an effective theory of massless degrees of freedom consisting of the 

unbroken subgroup associated with the r 3 = ( J J j generator of the SU(2) 

group. This model is then nothing else but the r 3 QED model of 0, discussed 
above, with the important difference that now the statistical unbroken gauge group 
U c (l) is compact. 

Compact groups are known to have non-perturbative monopole configurations 3 , 
which in (2+l)-dimensions are like instantons, having a point-like structure. Such 

3 Recently, in the slave-boson framework, compact QED3 theories have also attracted attention 
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monopole instantons are responsible in general for the generation of a small non- 
perturbative mass of the statistical photons of the U c (l) group [2T] : 

photon mass = m 7 ~ e~ 2So (2-12) 

where the above computations have been performed in a dilute instanton gas, Sq 
is the one-instanton action, and the factor of 2 in the exponent is due to the fact 
that the theory has fermions [2H]- The presence of a non-perturbative photon mass, 
implies that in the compact case there is no longer a massless pole in the electric 
current- current correlator (|2.5j) . As a consequence, the KT T = superconducting 
quantum-critical line of the non-compact case is absent, and the AF phase is sepa- 
rated by the (optimal doping) superconducting one in the phase diagram of fig. El 
by a pseudogap phase, which in this case extends all the way down to T = 0. 

3 Counting degrees of freedom: a tricky business 

We would now like to address the issue of dynamical generation of the fermion 
mass per se in all the above cases, (A) and (B). As we have discussed, such a mass 
generation is crucial for the existence of the novel quantum-critical 'anomalous' 
phases. To this end, in this section we shall briefly review first the non perturbative 
arguments of ^Uj, according to which dynamical mass generation can only occur in 
a theory with Nf four-component spinors, with Nf < N c = 3/2, showing that the 
previous large-iV treatment JB] overestimated N c . If this result were true in our 
case, this would mean that the above-described QED 3 model with two four- 

component spinors (spinons) constructed as in fig. Q would never have a quantum 
phase with broken % symmetry, since the spinons would be massless. However, as we 
shall discuss below, the arguments of jTU] may not go through in (2+l)-dimensions 
for a variety of reasons. At this stage we would like to mention, however, that even 
if such arguments were assumed to be valid in (2+l)-dimensional gauge theories, 
nevertheless their application to the r 3 -QED model [2 does not select a critical 
number for mass generation of the electrically charged fermions, leaving this task 
to a detailed SD analysis or other studies. Below we shall explain in detail why 
this is so, but we will also provide arguments supporting the point of view that the 
analysis of JU] are not applicable to the (2+l)-dimensional QED3 model either. 

We commence our discussion by first going over the new proposed constraint on 
strongly coupled field theories, of ^U] The constraint was based on a counting of 

from a different perspective, specifically in connection with the role of monopoles on the confining 
phase, and their implication on the Mott insulating phase of doped AF |22| . 
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massless degrees of freedom between IR and UV fixed points, leading to inequal- 
ities between the respective values of the renormalized free energy density. This 
resembles, but is not identical to, the celebrated C-theorem of Zamolodchikov for 
two-dimensional conformal field theories [24J. The conjectured constraint of [TU] . 
has not been proven; nevertheless, its validity has been verified in a number of phys- 
ically relevant examples. The conjecture of ^U] implies a reduction in the massless 
degrees of freedom as the theory flows under Renormalization Group (RG) from UV 
to IR fixed points. In what follows we shall therefore examine the applicability of 
the conditions leading to the constraint for the above- described QED3-like models 
of potential interest to the physics of doped antiferromagnets. 

For instructive purposes, let us first state the form of the constraint of ref. ^Uj 
and the main assumptions leading to it: The main thrust of the constraint is as- 
sociated with an inequality between the values of an extensive quantity, such as 
the free energy, at the UV and IR fixed points of asymptotically free field theories. 
Specifically, it has been argued that the quantity /uv = — hmT-^oo^-^ where T is 
the free energy of the system, and T is the temperature, playing here the role of 
a varying RG scale in the problem, counts the massless degrees of freedom in the 
UV, while a similar quantity: /ir = — lim T ^ ^^ counts the massless IR degrees 
of freedom of the system. The conjectured inequality states that: 

/m < /uv (3.1) 

An important ingredient for the validity of Eq. ()3.1|) is the asymptotic freedom of 
the theory in question. In fact, there are known examples that do not satisfy the 
inequality if they exhibit a non-trivial UV fixed point [10J. 

The inequality (|3.1|) has been applied to QED 3 , which is known to exhibit chiral 
symmetry breaking in the IR. The theory has a dimensionfull coupling, and in ^U] , 
it has been assumed free in the UV. According to JUL this implies the following 
contributions of massless degrees of freedom to /uv (i n three space time dimensions): 
one from the massless photon in (2+l)-dimensions, and (3/4) AN from N four- 
component free fermions. Hence one has at the UV: 

/ uv = 1 + 3N (3.2) 

On the other hand, in the IR, there is a dynamical generation of fermion masses, 
which implies a breaking of the global chiral symmetry of the massless theory 
U(2N) — > U(N) x U(N). In this case there are 2N 2 massless Goldstone bosons. 
The authors of [Till deal with non compact QED 3 , and hence the photon remains 
massless in the IR, which implies the following contribution to /ir : 

/ IR = 1 + 2N 2 (3.3) 
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From the inequality (|3.1j) . and (|3.2|) . (J3.3)) . then, one obtains a critical number 
of fermion flavours for chiral symmetry breaking to take place; in particular iV < 
N c = 3/2, a result, which, if true, would imply that the early papers on dynamical 
mass generation ^H] overestimated the critical number, placing it in the region 
3 < iV c < 4. 

Now let us go over these assumptions one by one: 

• (1) Asymptotic freedom: In the context of perturbation theory QED3 is super- 
renormalizable (i.e. it contains a finite number of divergent diagrams); this is 
a direct consequence of simple power-counting, given that the dimensionfull 
gauge coupling scales as [M] 1 / 2 . However, QED 3 is not asymptotically free, 
in the same way that the super-renormalizable (0 3 )4 is not (this latter theory 
has a dimensionfull coupling as well): the coupling does not go to zero for 
large momenta p — > 00. In fact, strictly speaking, the coupling does not 
display any energy dependence at all, because the corresponding (3 function 
vanishes, i.e. /3qed 3 = A* (de/d/j,) = 0. This is true simply because the vacuum 
polarization in d = 3 is finite, and thus no wave-function renormalization 
for the photon is needed. The above renormalization group equation simply 
states that e = e , where e is the initial tree-level coupling. Clearly such 
a theory cannot be asymptotically free (at least in a perturbative context), 
unless eo = 0; but then this would imply that the theory is non-interactive 
(i.e. trivial) throughout. This situation is to be contrasted to what happens 
in the case of a bona fide asymptotically free theory, such as QCD 4 . In this 
case /3qc*d 4 = —bog 3 (bo > 0); when solving this equation one obtains the usual 
logarithmic dependence of the coupling g on p, which, in the limit p — > 00, 
leads to an asymptotically non-interactive theory (i.e. the value of the coupling 
goes to zero, not just to some constant value). 

Given the discussion above, in order to obtain some sort of non-trivial structure 
in the UV, one must study QED 3 outside the realm of standard perturbation 
theory. In particular, arguments supporting asymptotic freedom of QED 3 , and 
hence the counting in the UV leading to (|3.2|) . are based on large N treatments. 
Indeed, in such an approximation one defines a dimensionless effective 
running number of flavours, which turns out to scale with the momenta as (in 
the non-local gauge): 

(effective coupling) 2 /scale = ~ x - ^ p/(y ; So - ^ (3.4) 

An obvious drawback of the large N treatment in this context is the fact that, 
even though the expansion is in principle systematic due to the presence of a 
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parameter (N) which may be formally taken to infinity, in practice the desired 
effects are obtained by stretching the validity of this approximation into a 
range of N values which one might consider as "uncomfortably small". In 
particular, an upper bound on the critical number of flavours N < N c ~ 3/2 
is obtained, beyond which chiral symmetry breaking is not possible. In this 
scheme the fundamental photon-fermion vertex is assumed to have the form 
T M = j fl A n (p), where A(p) is the wave function renormalization. For the case 
N = 2 there is no running and the model is not asymptotically free. For 
N > 2 the theory is asymptotically free, while for N < 2 the theory behaves 
like QED4 with a Landau pole. 4 



e(p) 








Figure 7: The dotted line represents the constant (non-running) coupling obtained 
from the /3qed 3 = condition, and the solid line is the running obtained through 
the effective charge of Eq. JRTTty . In both cases the theory is not asymptotically free, 
because the coupling saturates at a non-zero value. For comparison, the QCD^ cou- 
pling (perforated line) approaches asymptotically the value zero (asymptotic freedom) 

The Effective charge: 

Given that, due to the super-renormalizability of QED 3 one obtains no stan- 
dard running for the coupling, one might be tempted to explore other field- 

4 In Kondo and Nakatani (first item in the analysis was done in the Landau gauge, while 
in the Kondo-Murakami paper (second item in |25p an improved analysis was performed in the 
non-local gauge, where the vertex assumes the form T M = 7 /J G(p 2 , q 2 , k 2 ) but A{p) can be set 
identically equal to 1. 
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theoretic alternatives. In particular one may ask what sort of "running" one 
obtains from the "effective charge" in the case at hand. 

In QED4 the infinite subset of radiative corrections summed in the Dyson 
series generated by the one-particle-irreducible vacuum polarization ITr(p 2 ) 
defines an effective charge which is gauge-, scale-, and scheme-independent 
to all orders in perturbation theory: 

e 2 e 2 

elsip2) = 1 + iW) = TTW) ' (3 ' 5) 

where we have used e 2 R = (Z 2 Z 3 /Zf)e 2 and 1 + Ur = Z 3 (l + U) together 
with the QED Ward identity Z x = Z 2 to write e 2 ff (g 2 ) purely in terms of 
bare quantities. At p 2 /mj — > 00, the effective charge e 2 s (p 2 ) matches on to 
the running coupling e 2 {p 2 ) defined from the renormalization group: at the 
one-loop level, 

Km e 2 cS ( P 2 ) = eV) = g2 4 — - , (3.6) 

P 2 M-oo ! _ _££_ n/ l 0g (p2/ m 2) 

where rif is the number of fermion flavors. What makes the effective charge 
a particularly useful concept is its non-trivial dependence on masses, through 
the analyticity properties of H(p 2 )- In particular, when p 2 > 4m 2 , the vacuum 
polarization has imaginary part Qmll(p 2 ) given by 




SmlTV) = $(p 2 - 4m 2 )- 1 + — ± W 1 ^ , (3.7) 



By virtue of analyticity the real part KellQo 2 ) may be reproduced from Qmll(p 2 ) 
by means of a once-subtracted dispersion relation Thus, for the one-loop con- 
tribution of the fermion /, choosing the on-shell renormalization scheme, 

jtdV) = ay r ds ^i (3.8) 

Vtt/ J 4m 2 s{s-p 2 ) 

Finally, the one-loop SmllQo 2 ) is directly related, via the optical theorem, to 
the tree level cross sections for the physical processes e + e~ — > f + f~, with 
/ ± e, 

QmU(s) = a(e + e- -> /+/") (3.9) 

The aforementioned properties of the 4-d effective charge go through in the 
case of QED 3 , with the appropriate adjustments due to the fact that no renor- 
malization is needed (for example, the 3-d analogue of the dispersion relation 
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in Eq. (|3.8|) needs no subtraction) 5 . In the case of QED 3 the dimensionless 
effective charge reads 

elM,a= iT^m (3 - 10) 

where n(p 2 ) is the scalar co-factor of the one- loop vacuum polarization U^ u (p) = 
(p 2 9iiv ~ PfiPu)^-(p 2 ) in d — 3; in particular, H(p 2 ) = —a/8p. The qualitative 
behavior of e e g (p) is shown in fig.©, solid line: For large p the effective charge 
saturates to a non-vanishing value eo, i.e. it does not display asymptotic free- 
dom. 

• (2) IR infinities and the Mermin- Wagner theorem: 

The argument leading to the inequality ()3.1|) is based on a finite temperature 
extension. However, in two spatial dimensions at finite temperature there are 
IR infinities as T — * 0, which invalidate the well-defined nature of Goldstone 
bosons. Such IR infinities can be seen in the effective IR cut-off dependence of 
the running flavour number at any finite temperature N = iV(lnT) which stems 
from solving large- N SD equations in the presence of an IR cutoff [23 12H EH] • 
One may question the robustness of large- N treatments, but the presence of 
IR infinities is probably indisputable. 

The absence of Goldstone bosons, due to IR infinities (Mermin- Wagner the- 
orem), implies the following counting of massless degrees of freedom in the 
IR: fin = 1. Hence in that case the inequality (J3.1j) becomes formally (even 
assuming asymptotic freedom): 

fm = K /uv = 1 + 3iV -> N > 1/3 (3.11) 

hence we obtain no non-trivial information from this constraint. In fact the 
above counting is formal, given that the IR infinities imply InT divergences 
as T — > in the effective action, and hence in /tr. One therefore has to do 
the computation explicitly to regularize these divergences (201 |2H| , which leads 
to IR-cutoff- (and hence temperature-) dependent critical number of fermion 
flavours N C (T). 

• (3) This and the following items are physical reasons why the inequality (J3.1)) 
does not apply to systems of interest in our case. 

5 The concept of the effective charge has been generalized in a non-Abelian context through 
the use of the PT |261 177| . Recently it has been proposed that the effective charges provide a 
natural framework for the reliable study of the impact of threshold effects of the unification of 
gauge couplings in Particle Physics models [2H| 
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In most theories of condensed matter for high temperature superconductors 
there are four fermion interactions present, in addition to the gauge minimal 
coupling terms. It is well known that four-Fermi theories in less than 
four dimensions are renormalizable (in a large N treatment at least), and 
they exhibit a non-trivial UV fixed point. This is sufficient to invalidate the 
counting of degrees of freedom leading to /uv in this case, and hence the 
inequality (|3.1j) for such systems. 

• (4) Another important physical system which may imply a different critical 
number even if one accepts the inequality (|3.1|) is the model for high temper- 
ature superconductivity of [2 known as r 3 QED 3 , whose lagrangian is given 
in (J2.3|) . The important difference of this model from QED 3 (J2.7|) lies in the 
original global symmetries. Namely, in the r 3 QED model (|2.3|) . due to the op- 
posite couplings of the two fermionic 'colors' (each color is a four component 
spinor), the symmetry that is broken is a global fermion number symmetry 
J/j, = i'cjl^cj 6 - As discussed previously the symmetry is broken through 
the anomaly graph of figure El The figure denotes the S-matrix element 

< a^J^O >= sgn(m)e^ p p p / v /p^ (3.12) 

(where po is the energy) in the phase where there is mass generation m for 
the fermions. From the breaking of the global U(l) symmetry there is one 
Goldstone boson but again there is no local order parameter. In that case the 
counting of the inequality ()3.1|) implies 1 + 3iV > 2 i.e. N > 1/3. Note that 
in this model, if one applies the Mermin- Wagner theorem, according to which 
there is no well defined Goldstone boson, then the above inequality would only 
imply N > 0. At any rate, the coupling of the system to electro magnetism, 
promotes the global fermion number U(l) symmetry to a local electromag- 
netic symmetry, whose breaking results in superconductivity. In that case 
the would-be Goldstone boson is eaten by the longitudinal component of the 
electromagnetic potential, which now acquires a mass |2|. 

This completes our discussion on the limitations of the applicability of the con- 
straints of ref. PU] on the three-dimensional gauge systems of interest to us here. In 
the next section we proceed to review a SD analysis, which sheds some light on the 

6 It must be stressed at this point that one may define a chiral-like symmetry in this case 
generated by C75, where C is the (statistical) charge conjugation operator, however due to the 
presence of the discrete operator C this cannot be represented as a global U (2N) chiral symmetry, 
whose breaking leads to Goldstone bosons. At any rate, this symmetry is broken explicitly in 
the specific model by the coupling of the charged fermions (holons) to the real electromagnetic 
potential 0. 
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non-asymptotic freedom of the effective QED 3 model, the absence of a critical num- 
ber of fermion flavours for chiral symmetry breaking, but the existence of a region of 
the effective charge where the phenomenon takes place. In addition we also discuss 
the non-trivial IR fixed point structure of QED3. We believe that all these fea- 
tures are physically important; especially the last one implies the non- Fermi liquid 
behaviour of the relativistic spin liquid of nodal excitations under consideration. 



4 SD equations and Critical Number of Flavours 



In this section we will present a study of the issue of mass generation and chiral 
symmetry breaking in the framework of the SD equations. 



Figure 8: The Schwinger- Dyson (SD) equations for the photon and fermion self 
energies, and the vertex function. The blobs indicate the full (non-perturbative) 
corrections. 

The derivation of the SD equations for the photon propagator A^, the elec- 
tron propagator Sf, and the photon-electron vertex in QED3 proceeds following 
standard methods The result is schematically depicted in fig. |HJ 





+ 




(4.1) 
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where A ^(q) = q 2 g^ u + (f 1 - ^(p) = - m, and p x + p 2 + P3 = 0. 

The ellipses on the right-hand sides denote the infinite set of terms containing the 
two-particle irreducible four-point function [211 E2|- Although we are not working in 
the context of a large- N analysis, we note that the above truncation is compatible 
with working to leading order in resummed 1/N expansion. 
We next define the scalar quantities A, B and Q as follows: 

Sf ^ = A(k) f A ^ k "> = B(k)k 2 ' T ^ Pl > P2 > p ^ = GiPi'P^Psht* ( 4 - 2 ) 

Note that the form of Eq. (|4.2|) implies that the longitudinal pieces of the photon 
propagator will be discarded in what follows; this is motivated by the PT analysis 
presented in the Appendix, particularly points (iv) and (v) of the last subsection. 

4.1 The semi-amputated vertex 

Following |T2j and [331 we define the semi-amputated vertex G as 

G(pi,P2,Pa) = Z(jPi,P2,P3)G(pi,P2,Ps) ( 4 -3) 

with 

Z( Pl ,p 2 ,p 3 ) = B- l l 2 {p l )A- 1,2 {P2)A- l l 2 {p,) (4.4) 

This definition proves very useful in reducing the complexity of the set of SD equa- 
tions. In addition, as explained in |3IJ the quantity 



9R{Pi,P2,P3,)=eG(p 1 ,p2,P3) (4.5) 

provides a natural generalisation of the concept of the running or "effective" charge 
in the context of super-renormalizable gauge theories, such as QED 3 . 

The equation for the semi-amputated vertex G may be obtained from the third 
equation in (j4.1j) by multiplying both sides by the factor Z(pi,p 2 ,P3), i.e. 

/d 3 k - 1 1 1 

(27F GV ^ 7 ^ 7 "(fcT^F (4 - 6) 

where G 3 = G(ps, k + p 2 ,Pi — k)G(pi, —k, k — pi)G(p 2 , k, —k — p 2 ). Restricting 
ourselves to the case where the photon momentum is vanishingly small, one is left 
with a single momentum scale p. One can then define a renormalization-group f3 

function from this "running" coupling G(p) by setting (3 = p (dG(p) / dpj . In order to 

further simplify the SD equation for G(p) we make the additional approximation that 
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G 3 = G 3 (k), i.e. a cubic power of a single G(k) depending only on the integration 
variable k. 

Carrying out the 7-matrix algebra in d = 3- dimensional Euclidean space, one 
obtains: 

1 f d 3 k ~ 1 

G(p) ^ W + _ e2 |_ G 3 W __ (4 . 7) 

We observe that Z(p) — > 1 for p — > 00, where perturbation theory is valid. This 
is expected from the fact that in such a case the functions A(p), B{p) — > 1 trivially. 
Moreover, from ()4.7|) we see that, if G stays positive, which is expected for any 
physical theory, then, as a result of the positivity of the integrand, G(p) > 1 for 
any p. Thus, one has the following basic properties of G(p), which stem directly 
from the integral equation (j4.7|) : G{p) > 1, (for all p), and G{p) — > 1, p — ► 00. 
Assuming that in the IR regime, fc/a << 1, the inhomogeneous term Z(p) goes to 
zero, (this assumption has been justified by explicit calculations in [3*1] ) one can 
decouple the equation for the amputated vertex from the equations for A(p), B(p). 
Thus, one arrives at the homogeneous integral equation 

6 ^ = ¥lwf &(k) whw (4 ' 8) 

involving only one unknown function, namely G, which must be self-consistently 
determined. Note that Eq. (j4.8j) is invariant under the rescaling G — > G/e. This 
indicates a straightforward extension of the analysis to a large- N treatment, given 
that iV appears only as a multiplicative factor. Eq. (|4.8|) does not admit physically 
acceptable solutions, i.e. solutions with G > and finite. Indeed, setting p = one 
obtains after the (trivial) angular integration 

Finiteness of G(0) requires that the integrand of the right hand side of (|4.9|) converges 
at y — > and 00. The UV limit does not present a problem, because the kernel 
vanishes like y~~ 2 , which is consistent with the super-renormalizability of the theory 
as well as the fact that the amputated vertex tends to 1. In the IR limit y — > 
0, however, the kernel blows up. For the integral to remain finite at that point, 
as required by the finiteness assumption for G(0), G 3 (y) must approach zero as 
y a , a > 1/3, thereby implying that G(0) = 0. However for that to happen the 
integrand in ()4.9|) must change sign, which would in turn imply that G(y) itself 
must change sign somewhere in y. According to our assumption above this is not a 
physically acceptable situation. 



24 



The way to see that indeed the behaviour G(y) ~ y 1 ^ 2 as y — > 0, would be 
the only possibility is to convert the integral equation into a non-linear differential 
equation. To this end, we perform the angular integration in (|4.8|) . to arrive at the 
equation: 

2 a f°° dk A3 „^_ i k + p i 



where we have set e 2 = 8a to make contact with the usual large- N definition [TH] . 
For us, however, the number of fermion flavours is not assumed to be necessarily 
large. Introducing the dimensionless variables x = p/at and y = k/a, one obtains in 
the limit x << 1 

3vr 2 x y Q y y - x 3vr 2 V.a; 2 J ^ y 2 

(4.11) 

Differentiating appropriately with respect to x, and conveniently rescaling G by 
setting G = J j^tG, we arrive at the following differential equation for small x: 

x 3±}± + 3x ^ + G ^ x ) = 0, x «i (4.12) 
ax 2 ax 

The obvious solution of this equation, for x « 1, is G{x) = ^x 1//2 . Note that this 
solution would imply a 'trivial IR fixed point structure' given that its associated f3 
function vanishes at x — 0. As we shall demonstrate below this is in fact the only 
solution with G > 0. Indeed, upon the change of variables = 2x 2 , G = 2 3 / 4 r/(£), 
the equation becomes of Emden-Fowler type jSHHSSj: 



'2 



- r 3/ V(0 = o, e - +oo (4.i3) 

As discussed in the mathematical literature [HE], the only positive solution of 1)4 
as £ — > +oo, has the form as £ — > +oo: 

which is exactly the solution G(x) = ^x 1 / 2 found above, in the region x — > 0. 

The above analysis suggests that no non-trivial IR-fixed point is possible in QED 3 
in the absence of an IR cut-off, as already conjectured in ref. 
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4.2 The IR fixed point 

We next study the behaviour of the SD equation for G(p) in the presence of an IR 
cut-off. W shall consider the case of a fermion mass gap m(p) = Y,(p)/A(p), where 
E(p) is the fermion self energy. In that case the fermion propagator Sp 1 becomes: 



S F (k) 



A(k)Qf + m f (k)) 



(4.15) 



and we assume that rrif(p) ~ wi/(0) = rrif 7^ 0. In that case the integral equation 
()4.8|) becomes: 



G{p) 



e 
3" 



d 3 k 



G 3 (k) 



(2vr) 3 (k 2 + m 2 ){k - p)< 



+2mi 



d 3 k 



G 3 (k) 



(2it) 3 (k 2 + m 2 ) 2 {k - pf 



(4.16) 



Performing the angular integrations one arrives at: 



G{x) 



3n 2 x 



y + x 



y-x 



G\y) 



(4.17) 



where x = p/a, m = rrif/a are dimensionless, and 



(y 2 + m 



2\2 



(4.18) 



Differentiation with respect to x yields: 

2 



x— G(x) 
ax 



3n 2 x 



dyf(y) In 



y + x 



y-x 



+ 



2xy 
x 2 — y 2 



G 3 (y) 



(4.19) 



One observes that formally as x —>■ the right-hand-side vanishes, provided that 
G is finite. This indicates the existence of a fixed point. As we shall show below 
this is confirmed analytically by converting the integral equation into a non-linear 
differential equation. To accomplish this, one expands the logarithms for small 
x « 1, and then differentiates with respect to x, arriving at 



d 2 



d 



x(x 2 + m 2 ) 2 —G(x)+3(x 2 + m 2 ) 2 ^G(x) , , , 
dx A dx 67T Z 



x 



3m 2 )G 3 (x) = (4.20) 



In the IR region x « m the equation ()4.20j) is approximated by: 

-G 3 (x) = 



d 2 ~ d - 

x— G(x) + 3— G(x) + 



dx 2 



dx 



n 2 m 2 



(4.21) 
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and is immediate to see that a special power-law solution is given by (for positive 
G(x))by: 

G(x) = m7 c-^=x- 1/2 (4.22) 
4a/2 

Notice the IR divergence of this type of solutions even in the presence of a (bare) 
fermion mass. The associated renormalization-group j3 function for this case reads: 

P(x) = -~G ~ x~ 1/2 -> +00, as x -> (4.23) 

indicating the absence of an IR fixed point. The associated operator appears to 
be relevant (negative scaling dimension), which implies the possibility of the theory 
driven to a non-trivial fixed point. 

However, in the IR regime x « 1, one can find a different type of solution: 



V3_ 

2y/2l + cV 



G = mir - V / _ 1 | ^ 2 , x -> (4.24) 



where c is a constant of integration to be fixed by the boundary condition at x = 
implied by the integral equation, to be discussed later on. For physical solutions c 
is assumed positive. This type of solutions has a renormalization-group /3-function 
of the form: 

(3 = -G(x) + J^-G 2 (x) ~ - X - 0, x - (4.25) 
Vo7rmc (1 + (rxj 

from which we observe the existence of a non-trivial (non-perturbative) IR fixed 
point at G* = 7r ™J^ c > 0. Such a fixed point is the result of the dynamical generation 
of a parity- invariant, chiral-symmetry breaking fermion mass |18J, indicating the 
connection of the phenomenon of chiral symmetry breaking in QED3 with a non- 
trivial IR fixed point structure. 

The non-trivial fixed-point solution (J4.24|) is compatible with the integral equa- 
tion (J4.16|) for some values of the fermion mass m to be specified below. Indeed, 
one can derive a boundary condition for G(0) from (j4.16|) . which reads: 

^4r^^p»j* ( 4, 6 ) 

In contrast to the massless case (j4.9j) . G(0) is now a finite constant, nmc a/3/8, 
as seen from 1)4.24)1 . and this allows for a compatibility of the solution ()4.24jl with 
([4.26)1 . provided that mG(0) satisfies certain conditions. Such conditions have been 
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derived in |34| . where a full analysis of the non linear coupled system of vertex, 
fermion and photon self-energy integral equations has been performed. The result 
is summarized in the equation: 

fit* . 



m < -^z: — tSh> > —r~ - 5 - 88 • ( 4 - 27 ) 

IT 



G(0)-^' 



This condition will act as a boundary condition for the allowed values of m in a 
mass-coupling diagram, as we shall discuss later on. 

4.3 Dynamical generation of the fermion mass gap 

In the previous subsection we have assumed the presence of a finite fermion mass, 
which we have treated effectively as an arbitrary parameter of the model. In this 
subsection we turn to the full problem, and study the dynamical generation of this 
mass, by deriving it self-consistently from the corresponding SD mass-gap equation. 
The equation for the gap reads 

/d 3 k 1 1 

m^Wm (428) 

where M(k) = H(k)/A(k) is the mass function, and we have pulled out factors of 
A{p) appropriately so as to be able to define an amputated vertex function Q(k). 
After standard algebraic manipulations, using dimensionless variables, in units of 
a = e 2 /8, M = M(k)/a, x = p/a, y = k/a, and working in the regime of low 
momenta x « 1, one arrives at the following differential equation: 

In the relevant region x 2 « M 2 « 1, we neglect x 2 next to M 2 in (|4.29J) and use 
the solution ()4.24|) for G{x) ~ M\/^ttc as x — » 0. The result is: 



d 2 ~ d ~ 

x—M(x) + 3— Mix) + 9c 2 M(x) =0, x -> (4.30) 
dx 2 dx 

from which one obtains a power series expression for the dynamical mass: 

00 i q 

M{x) = C,x- 1 Y,i~ l T ( 3c ) 2n+2 xl+n rrTTZ S ~ + °( x )> x ~+ 

nil (3 + n) 2 

(4.31) 
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From this one obtains the following relation between G(0) and M(0) = rrif/a: 

M (0) = m f /a = m==(—j Cl /3 &' 3 {0) (4.32) 

4.4 Comparison with the large N 

It is important to compare the above results with those obtained within the context 
of a large- A analysis [T%] . In particular, at first sight it seems that the relation 
(14.32)1 does not have a critical coupling, above which dynamical mass generation 
occurs. However, because the result 1)4.32)1 has been derived in the context of the 
solution (|4.24|) . one should bear in mind the restrictions characterizing this situation. 
In particular, as has been explained in detail in [33] the following conditions must 
hold: 

£5/ 3(0) _ 1^g 2 / 3 (0) (l - -) < 0, 6(0) > ^ (4.33) 

5 3 (12ttCi) 1/3 I V 5 

This restriction implies a critical coupling, G c = 12V6/5 ~ 5.88 but it is derived 
in a way independent of any large- A analysis. The way to understand ()4.33)) is the 
following: one should first fix a range of G(0), with G(0) > 5.88, and then use a 
C\ that will be such that, within this range of the couplings, eq. ()4.33)1 is satisfied 
for masses fa « 1. As can be readily seen, the bound for C\ obtained from the 
requirement that m << 1 is far less restrictive than the one associated with (j4.33j) . 
provided G(0) is not too close to the critical G c (0), where the mass m vanishes of 
course. For instance, for G(0) = 0(8), the upper bound on C\ from ()4.33j) is of 
order (9(10~ 4 ), while for G(0) = 6 the upper bound is C\ < 4. Notice that the 
bound is very sensitive to small changes in G(0). 

A typical situation is depicted in fig. El for two values of G\ = 1CT 5 , 10 -2 . We 
observe that the case C\ = 1CT 2 yields an upper bound in the mass which is of order 
0.8 and hence should be discarded on the basis that it is not small enough. On the 
other hand the value C\ = yields an acceptable upper bound m ~ 0.1. In that 
case, from fig. El we observe that the allowed region of m is 

0.08 <m< 0.12 (4.34) 

The corresponding regime of the couplings G(0) is: 

5.88 < 6(0) < 11 (4.35) 

The above results are to be compared with the corresponding regimes of masses and 
couplings derived in a large- A" analysis. We recall that, in the context of a large A 
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Figure 9: Fermion mass versus the IR-value of the coupling G(0) using \4-32j) 
(dashed curves), for two values of C\ = 1CT 5 fl ower dashed curve) and C\ — 10 -2 
(upper dashed curve). The continuous curve is viewed as a boundary condi- 

tion. The value C\ = 10~ 2 should be excluded on grounds of yielding too high values 
of the mass fh. 



treatment, and to leading order in 1/N resummation, the following solution for the 
dynamically- generated m is found [18J 



m ~ C(l)exp 




(4.36) 



where g 2 c = |j is the critical coupling, above which dynamical mass generation 
occurs PH]- I n fact g c is interpreted as an inverse critical number 1/N C of four 
component fermion flavours. When 1/N 2 or higher-order corrections are included, 
the dynamical mass generation procedure yields 3 < iV c < 4 (a result that as we have 
seen has been refuted by the inequality of [TO] , discussed in the previous section). 

Compatibility of the solution (j4.36J) with the constraints obtained from the inte- 
gral equations for the vertex implies jHl] the existence of an upper bound on fermion 
masses, m < m max , where m max is defined through the intersection of appropriate 
curves coming from the non linear constraints (c.f. fig. ITTHl . This yields m max ~ 0.3. 

On the other hand, for large momenta, we know that G — > 1. Physically one 
expects a monotonic decrease of G(x) over the entire range of x G [0,oo). This 
would occur in our case if and only if G(0) > 1, which, in the context of the large- N 
result of |18| I2]. implies a minimum bound for the fermion masses m > m min ~ 0.03. 
Actually, as we shall argue in the next section, G(0) should be comfortably larger 
than a/3/2 for self-consistency of our approximations. 
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Figure 10: Fermion mass versus the IR-value of the coupling G(0). The solid curve 
represents the condition derived from the integral equation for the vertex, whereas 
the dashed line represents the solution obtained from the standard gap equation in 
the large-N treatment. 



Hence, we see that the monotonicity of the running coupling can be achieved in 
the context of a large- N treatment, if the mass m lies in the following regime: 

0.03 < m < 0.3 (4.37) 

or equivalently if the coupling at the IR point G(0) is restricted in the regime: 

1 < G(0) < 2.5, (4.38) 

which is to be compared with (|4.35|) . 

From the physical point of view of applicability of the QED 3 models to the 
theory of high temperature superconductors, the above results imply constraints on 
the microscopic parameters of the underlying lattice models, whose long-wavelength 
limit is the three dimensional gauge theory under consideration. It must be noted 
in this case that the above analysis applies equally well to both QED 3 and r 3 -QED 
models. For instance, in the models of ^312] the effective (gauge-invariant) coupling 
G(0) may be expressed in terms of the parameters of the microscopic condensed- 
matter lattice systems, as: 

G(0) 2 ~4a- 5 ) (4.39) 

e z 

where S expresses the concentration of impurities in the system (doping), and J 
denotes the Heisenberg (antiferromagnetic) exchange energy. Hence, on account 
of ItOHjl . Eq. (jOHj) implies that 6 < (J/e 2 )(l -5) < 11 for superconductivity 
to occur. In phenomenologically acceptable models |2 ; e 2 /J ~ 0.1, which implies 
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an upper bound on 5 ~ 0.4, which is physically consistent. By slightly modifying 
this ratio we may even obtain also lower limits for 5, and hence a region where 
the phenomenon of fermion mass generation of the nodal excitations occurs. How- 
ever, the reader should bear in mind that the above- described limiting values are 
rather indicative at present, given that a complete quantitative understanding of 
the underlying dynamics of high-temperature superconductivity from an effective 
gauge-theory point of view is still lacking. 

5 Conclusions and Outlook 

In this review we have described some unconventional phases of QED 3 and related 
models, occurring at T = 0, including KT superconductivity. Such features of QED 3 
are important when discussing the phase diagrams of underdoped AF. For the novel 
phases to occur it is necessary to have a dynamical generation of fermion mass for 
the case of two four-component spinors, which corresponds to the physical case. 

This last feature has been questioned recently in JU], by resorting to a con- 
jectured inequality. Should this inequality be applicable in this context, it would 
strongly disfavour the generation of chiral symmetry breaking fermion mass for the 
case two four component spinors, contrary to the findings of various analysis based 
on SD equations. However, we have presented a number of arguments as to why, 
in our opinion, the conditions necessary for the applicability of the inequality in 
question are not fulfilled in the relevant models. 

One might expect that a future resolution of some of the above issues can be 
furnished by lattice models: QED3 with an even number of fermion flavours can 
indeed be simulated on the lattice [37], as it does not suffer from a fermion deter- 
minant sign problem. However, this is not the case for the non-compact r 3 -QED 7 , 
whose simulations on a lattice still remain a big challenge, due to a sign problem in 
the respective fermion determinant, as a result of the T3 structure of the statistical 
gauge fields. Nevertheless, the compact case of T3 QED, which can be considered as 
a broken phase of a SU(2) gauge theory, can be simulated. 

However, as remarked, in [3*0] . caution should be exercised when one simulates 
QED 3 models on the lattice, due to the IR cutoff-dependence of many of the quanti- 
ties entering the pertinent chiral symmetry breaking dynamics, and in particular the 
critical number of fermion flavours. This is because of the non-trivial IR divergences 
of the model, which need regularization by means of the scale /i 8 . 

7 We thank S. Hands for an informative discussion on this point. 

8 It is worth remarking at this stage that the dependence of N c on the IR cutoff mass scale, 
fi, is actually similar to its temperature dependence, discussed in [221, given that, under certain 
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In jHUj both a numerical study and an analytic (approximate) expression have 
been given for the dependence of N c on \i: (z/ c /2)ln(a//z)+2arctan(z/ c ) = ir, with v c = 
\J 32 / (iV c 7r 2 ) — 1, and a = Ne 2 /8, and e the QED 3 bare coupling (with dimensions 
of square root of mass). ^From either the above analytic formula or the numerical 
analysis presented in [3U] follows that, in order to find chiral symmetry breaking for 
N = 2, one should use at least a ratio ire 2 / fi = 5 x 10 3 . 

Note that in order to extract safe conclusions for chiral symmetry breaking on 
a lattice, with a dimensionless lattice coupling f3 ~ l/e 2 a, where a is the lattice 
spacing, it is imperative that the volume of the lattice L 3 should be large compared 
to any other dynamically generated correlations of the system. Thus, the physical 
volume (L/P) 3 must be large. If we set La ~ 7r//x [30j, then L/ (3 ~ 7re 2 //i, and hence 
to observe a chiral symmetry breaking for the case of two four- component fermion 
flavours (N = 2) one needs L//3 > 10 3 , according to the above remarks. However, 
typical lattice simulations in the existing literature [37] use L/(3 = (9(10 2 ). This, 
in turn, may explain why chiral symmetry breaking is not seen in such cases (for 
N = 2). 

Given the practical difficulties in reaching such high values of L/f3 = 0(5 x 10 3 ) 
(with the presently available computing power), one may want to explore the possi- 
bility of extracting exact results on some of the physical issues regarding the novel 
quantum critical phases discussed in the present article. One such possibility would 
be to embed the gauge theories in question into some super ■symmetric (SUSY) the- 
ory, with at least Af = 2 conserved SUSY charges, which is the minimum number 
of supersymmetries to allow for exact results in (2+l)-dimensions. This has been 
argued to be possible in the context of semi-realistic condensed matter situations in 
[To] , from the point of view of a composite operator field theory, based on the dynam- 
ics of composites consisting of spinons and holons. The basic idea is that at certain 
regions of the parameter space of appropriate microscopic condensed matter mod- 
els (such as extended t — j models with next-to- nearest neighbor interactions [35]). 
there is an effective supersymmetry between spinons and holons, characterising the 
dynamics of the effective continuum field theory near nodal points. This symmetry, 
originally an Af = 1 SUSY, can be elevated to an extended Af = 2 SUSY, as a re- 
sult of the low-dimensionality of the model, implying the existence of a topologically 
conserved current, which supplies the latter supersymmetry structure 

The reason for working with composites, instead directly with spinon and holon 
constituents, lies on the fact that in this way one obtains dynamically massless gauge 
fields (Abelian), as particular combinations of spinons and holons. The resulting 
Af = 2 SUSY transformations of the composite fields have been ensured up to quartic 

circumstances, one may associate the temperature with an IR cutoff. 
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order in the spinon and holon constituents in ^B], where we refer the interested 
reader for details. In this approach only the spinon (magnon) bosonic fields, z, and 
the fermionic holon fields, tp, are fundamental degrees of freedom. The statistical 
gauge fields are induced as particular combinations of these fields, consistent with a 
composite Af = 2 SUSY at specific regions of the parameter space of the microscopic 
condensed-matter model. 

The resulting effective continuum model turns out to be the Af = 2 SUSY Abelian 
Higgs model. There are some exact results associated with the phase structure of 
the model, in particular the topology of the so-called moduli space, i.e. the space 
of the (complex) parameters of the SUSY model. These results are detailed in 
[To] and we shall not discuss them here. However, we do note that they include 
a passa ge from a pseudogap to an unconventional KT superconducting phase for 
the composite model, in the compact gauge field case, of similar nature to the one 
discussed here for the constituent excitations of spinon and holons. This passage 
may be understood qualitatively as follows: the Af = 2 SUSY Abelian Higgs model 
(as is otherwise called the Af = 2 SUSY QED 3 (SQED)) possesses a Higgs phase, 
where the gauge field is massive; this phase has been argued to correspond to a 
pseudogap phase, which, in the compact SQED case, may also be characterised by 
stripe-like configurations, due to domain walls of the compact gauge field [T£]. The 
compact SQED contains monopoles with both negative and positive charges, and 
the corresponding antimonopoles. The fugacities of the monopole configurations of 
the compact case may depend on doping concentration in such a way that at certain 
doping concentrations the fugacities of the negative charge monopoles vanish, leaving 
only monopoles of positive charge (say +1), and antimonopoles of charge -1. Such a 
case is similar to the case where the compact U(l) gauge theory is embedded into an 
SU(2) gauge theory. In such SU(2)-like theories it can be shown [TH] that there are 
no stripe phases, and moreover the statistical photon remains exactly massless (22] , 
thereby implying the onset of a KT superconducting phase, as in |2j, and the end 
of a striped pseudogapped phase. 

Finally, in [3D] the issue of dynamical mass generation in Af = 1 SUSY QED3 
was studied using superfield SD equations. It was shown that the presence of 
a supersymmetry-preserving mass for the matter multiplet stabilizes the infrared 
gauge couplings against oscillations present in the massless case, thus inferring that 
the massive vacuum is dynamically selected at the level of the quantum effective 
action. 

It must be noted that such supersymemtric points in general do not correspond 
to realistic values of the parameters of the underlying condensed-matter microscopic 
model. Nevertheless, they may be viewed as implying the possibility of using such 
extended supersymmetry (between spinon and holon composites) as a tool for ex- 
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tracting exact analytic information on the phase diagram of doped antiferromagnets 
as follows: one studies first such extended SUSY points, and then breaks the M = 2 
SUSY explicitly down to M = by varying the parameters of the model in order 
to approach the physical regime. One may then hope to construct models in which 
such a breaking would result in a situation where the SUSY partners of the physical 
excitations acquire masses higher than the highest mass scale in the problem (e.g. 
Debye screening), and hence can be safely discarded. Such issues definitely deserve 
further studies, and indeed may pave the way for obtaining an exact analytic un- 
derstanding of the complicated phase diagram of doped antiferromagnets, without 
resorting to lattice simulations of the models. It remains to be seen whether such 
hopes can be realised within the context of phenomenologically realistic condensed- 
matter theories. 
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Appendix: The Pinch Technique 

It is well known that off-shell Green's functions depend in general on the gauge- 
fixing procedure used to quantize the theory, and in particular on the gauge-fixing 
parameter (GFP) chosen within a given scheme. The fermion self-energy for 
example, is GFP-dependent already at the one-loop level. The dependence on the 
GFP is in general non-trivial and affects the properties of a given Green's function. 
In QED 4 , and in the framework of the covariant gauges, depending on the choice of 
the GFP £, one may eliminate the UV divergence of the one-loop electron propagator 
S(p, £) by choosing the Landau gauge £ = 0, or the IR divergence appearing after 
on-shell renormalization by choosing the Yennie- Fried gauge £ = 3. The situation 
becomes even more complicated in the case of non-Abelian gauge theories, where 
all Green's functions depend on the GFP. Of course, when forming observables the 
gauge dependences of the Green's functions cancel among each other order by order 
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in perturbation theory, due to powerful field-theoretical properties, a fact which re- 
duces their seriousness. However, these dependences pose a major difficulty when 
one attempts to extract physically meaningful information from individual Green's 
functions. This is the case when studying the SD equations; this infinite system of 
coupled non-linear integral equations for all Green's functions of the theory is inher- 
ently non-perturbative and can accommodate phenomena such as chiral symmetry 
breaking and dynamical mass generation. In practice one is severely limited in their 
use, and a self-consistent truncation scheme is needed. The main problem in this 
context is that the SD equations are built out of gauge- dependent Green's functions; 
since the cancellation mechanism is very subtle, involving a delicate conspiracy of 
terms from all orders, a casual truncation often gives rise to gauge- dependent ap- 
proximations for ostensibly gauge-independent quantities (HH HO] • The study of SD 
equations, and especially of "gap equations" , has been particularly popular in QED 
jH], and even more so in QCD 02], where it has been intimately associated with 
the mechanism that breaks the chiral symmetry. Similar equations are relevant in 
QED 3 , where the IR regime of the theory is probed for a non-trivial fixed point [2], 
for technicolor models [13], gauged Nambu-Jona-Lasinio models |H], and more re- 
cently color superconductivity |%5| . A similar quest takes place in top-color models, 
where the mass of the top quark is generated through a gap equation involving a 
strongly interacting massive gauge field jlH] • The usual conceptual drawback is that 
sooner or later one is forced to choose a gauge, resorting to a variety of arguments; 
but gauge choices cast in general doubts on the robustness of the conclusions thusly 
reached. 

To address the problems of the gauge- dependence of off-shell Green's functions 
a method known as the PT has been introduced fTJ ^] . The PT is a diagrammatic 
method which exploits the underlying symmetries encoded in a physical amplitude 
such as an ^-matrix element, or a Wilson loop, in order to construct effective Green's 
functions with special properties. The aforementioned symmetries, even though they 
are always present, they are usually concealed by the gauge-fixing procedure. The 
PT makes them manifest by means of a fixed algorithm, which does not depend on 
the gauge-fixing scheme one uses in order to quantize the theory, i.e. regardless of 
the set of Feynman rules used when writing down the S'-matrix element. The method 
exploits the elementary Ward identities triggered by the longitudinal momenta ap- 
pearing inside Feynman diagrams in order to enforce massive cancellations. The 
realization of these cancellations mixes non-trivially contributions stemming from 
diagrams of different kinematic nature (propagators, vertices, boxes). Thus, a given 
physical amplitude is reorganized into sub-amplitudes, which have the same kine- 
matic properties as conventional n-point functions and, in addition, are endowed 
with desirable physical properties, such as GFP-independence. Finally, the PT 
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amounts to a non-trivial reorganization of the perturbative expansion The role of the 
PT when dealing with SD equations is to (eventually) trade the conventional SD 
series for another, written in terms of the new, gauge-independent building blocks 
[TTJ EH H7j The upshot of this program would then be to truncate this new series, 
by keeping only a few terms in a "dressed-loop" expansion, and maintain exact 
gauge-invariance, while at the same time accommodating non-perturbative effects. 
We hasten to emphasize that the aforementioned program has not been completed; 
however, a great deal of important insight on the precise GFP-cancellation mecha- 
nism has been accumulated, and the field-theoretic properties of gauge-independent 
Green's functions have been established in detail. 

An explicit one-loop example. 

We next explain how the PT gives rise to effective, gauge-independent fermion 
self-energies at one- loop, for the case of QED and QCD |48j . As will become clear 
in what follows, the procedure described does not depend on the dimensionality of 
space-time; in particular, it applies unaltered at d — 3, 4. We will assume that 
the theory has been gauge-fixed by introducing in the gauge-invariant Lagrangian a 
gauge-fixing term of the form ^(d^A^) 2 , i. e. a linear, covariant gauge; the parameter 
£ is the GFP. This gauge-fixing term gives rise to a bare gauge-boson propagator of 
the form 

r £ £ 



2 



(5.1; 



which explicitly depends on £. The trivial color factor 5 a t> appearing in the (gluon) 
propagator has been suppressed. The form of A A4I/ (€, £) for the special choice £ = 1 
(Feynman gauge) will be of central importance in what follows; we will denote it by 
AJ„W, i.e. 

A, u (£,l)=A^(£) = -^ gfll/ . (5.2) 
A^ and A^{£) will be denoted graphically as follows: 

ttooooott = zA Mi ,(£, £), -www = iAf^ u (£). 

For the diagrammatic proofs that will follow, in addition to the propagators A^ u {£) 
and A^(£) introduced above, we will need six auxiliary propagator-like structures, 
as shown below: 
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•Jl_ 
1 1 



•wxfwv = -j± 



<WVWW= -zt #wvw\» = 4r 

t t 

All of these six structures will arise from algebraic manipulations of the original 
A^ v {€). For example, in terms of the above notation we have the following simple 
relation (we will set A = £ — 1): 

TTOOOOOC = WWW -| 

We next turn to the study of the gauge- dependence of the fermion self-energy (elec- 
tron in QED, quarks in QCD). The inverse electron propagator of order n in the 
perturbative expansion has the form (again suppressing color) 

S£ 1 (p,fl=*-™-£ w (p,fl (5-3) 

where E^(p, £) is the n—th order self-energy. Clearly £(°) = 0, and Sq 1 (p) = tf — m. 
The quantity H^ n >{p, £) depends explicitly on £ already for n = 1. In particular 

£ {1) (M) = | [^/Sofr + Q/VtU) = S«(p) + AE«(p) (5.4) 

£%) = EW(p, 1) = J [d£]YS (p + £)YK^) ( 5 - 5 ) 
and 

= j [ ^S {p + t)S Q \p)-S,\p) (5.6) 

In the above formulas [cff] = g 2 /i 2e , with D = 4 — 2e the dimension of space- 
time, /i the 't Hooft mass, and g the gauge coupling (g = e for QED, and g = g s 
for QCD). The subscripts "F" and "L" stand for "Feynman" and "Longitudinal", 
respectively. Notice that £^ is proportional to Sq 1 (p) and thus vanishes "on-shell". 
The most direct way to arrive at the results of Eq. (|5.fi|) is to employ the fundamental 
WI 

f = S \p + £)-S \p), (5.7) 



with 
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Figure 11: One loop diagram contributing to the QED/QCD fermion self-energy. 

which is triggered every time the longitudinal momenta of A M „(£, £) gets contracted 
with the appropriate 7 matrix appearing in the vertices. Diagrammatically, this 
elementary WI gets translated to 




Then, the diagrammatic representation of Eq. (j5.4j) . Eq. (|5.5j) . and Eq. (j5.6p will be 
given by 




(5.8) 

When considering physical amplitudes, the characteristic structure of the longi- 
tudinal parts established above allows for their cancellation against identical con- 
tributions originating from diagrams which are kinematically different from fermion 
self-energies, such as vertex-graphs or boxes, without the need for integration over 
the internal virtual momenta. This last property is important because in this way 
the original kinematical identity is guaranteed to be maintained; instead, loop inte- 
grations generally mix the various kinematics. Diagrammatically, the action of the 
WI is very distinct: it always gives rise to unphysical effective vertices, i.e. vertices 
which do not appear in the original Lagrangian; all such vertices cancel in the full, 
gauge-invariant amplitude. 
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To actually pursue these special cancellations explicitly one may choose among 
a variety of gauge invariant quantities. For example, one may consider the current 
correlation function defined as (in momentum space) 

hM = * J d 4 xe^'(0|T [J M (ar) J v (0)] |0> = (g^q 2 - q^)I(q 2 ) , (5.9) 

where the current J^{x) is given by J^{x) = :Q(x) r y fJ ,Q(x) :. I^ v {q) coincides with 
the photon vacuum polarization of QED. 

To see explicitly the mechanism enforcing these cancellations in the QED and 
QCD cases, we first consider the one-loop photonic or gluonic corrections, respec- 
tively, to the quantity I^ u . Clearly either set of corrections is GFP-independent, 
since the current J^{x) is invariant under both the £7(1) and the SU(3) gauge trans- 
formations. 

The relevant diagrams are those shown in FiglTTl To see the appearance of the 
unphysical vertices, we carry out the manipulations presented in Eq. (j5.4j) . Eq. (j5.5|) . 
and Eq. (j5.6|) . or, equivalently, in Eq. (j5.8|) . this time embedded inside I^q). Thus, 
from diagrams (b) and (c) we arrive at 




We thus see that since the action of the elementary WI of Eq. (|5.7jl amounts to 
the cancellation of internal propagators, its diagrammatic consequence is that of 
introducing an unphysical effective vertex, describing an interaction of the form 
llQQ or iGQQi depending on whether we consider photonic or gluonic corrections. 
This type of vertex may be depicted by means of a Feynman rule of the form 




being \i the index of the external current. 

To see how the above unphysical contributions cancel inside we turn to 
diagram (a). The action of the WI may be translated to the following diagrammatic 
picture 
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It is then elementary to establish that the two diagrams on the right-hand side of 
the above diagrammatic equation add up. 

Summing up the two equations above, it is clear how the gauge dependent part of 
the one loop amplitude cancel completely. Having proved that the GFP-dependent 
contributions coming from the original graphs containing E^^(p, £), i.e. Fig JTTT b) 
and Fig lTTT c) cancel exactly against equal but opposite propagator-like contribu- 
tions coming from Fig lTiT a). one is left with the "pure" GFP-independent one-loop 
fermion self-energy, E^(p). Clearly, it coincides with the E^P(p) of Eq. (|5.5|) . i.e. 

sW(p) = 4 1) (p). (5.10) 

All-order results. 

The generalization of the PT to all orders has been recently accomplished |4*^ IBTI]. 
The main points are the following: 

(i) In a fully non-Abelian context the longitudinal momenta responsible for the 
various rearrangements between different Green's functions do not stem solely from 
the tree-level propagator of the gauge boson (as in the previous one-loop example) 
but also from the vertices carrying momenta, i.e. the three-boson vertex Tq^ ' [T2] . 
It turns out that this latter type of longitudinal momenta triggers a fundamental 
cancellation taking place between graphs of distinct kinematic nature (s-channel 
versus t-channel graphs), shown to lowest order in In particular, when the s- 
channel and t-channel diagrams of FigIT2j) are contracted by a common longitudinal 
momentum, one obtains from either graph a common, propagator-like part, which 
eventually cancels against the other 26J. These parts display the characteristic 
feature that, when depicted by means of Feynman diagrams, they contain unphysical 
vertices (Figll2|). i.e., vertices which do not exist in the original Lagrangian, and 
cancel in any observable quantity. 

(ii) The generalization of the above cancellation to all orders in perturbation 
theory has been demonstrated recently in |4TH 130] . In the all-order case one con- 
siders the divergence of the four-point function A a ^A b u q l q> , with the gluons A a ^ A b v 
off-shell, and the quarks q 1 , (p on-shell. The aforementioned four-point function 
constitutes a common kernel to all self-energy and vertex diagrams appearing in 
the process q m q a — > q l qi (note that the diagrams of Fig^] are simply the tree- level 
contribution to the amplitude A^ A h v q % qi) . As has been shown in OED1 the judi- 
cious exploitation of the all-order Slavnov- Taylor identity that this Green's function 
satisfies allows for the all-order generalization of the PT procedure. 

(iii) The effective gauge- independent Green's function constructed by the PT 
coincide to all orders with the the background field method Green's functions when 
the latter are computed in the (quantum) Feynman gauge [HU ED] ■ 
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Figure 12: The tree-level version of the fundamental s-t channel cancellation. 

(iv) The equality of Eq (|5.1()jl persists to all orders: after all gauge-cancellation 
have been carried out, the gauge- independent PT fermion self-energy coincides with 
the conventional one computed in the renormalizable Feynman gauge. The gen- 
eral statement of (iii) is of course valid in this well, simply because the 
fermion self-energy computed in the renormalizable Feynman gauge happens to be 
identical (to all orders) to the fermion self-energy computed in the background field 
method Feynman gauge. This last statement is of course not true in general; Green's 
functions computed in the renormalizable Feynman gauge do not coincide with the 
corresponding Green's functions computed in the background field method Feynman 
gauge. What is always true however is (iii). 

(v) In QED (but not in QCD) the statement of (iv) is true also for the off-shell 
photon-fermion vertex T^. 
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